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The human influenza A (H3N2) virus dominated the 2014-2015 winter season in many countries and caused massive morbid- 
ity and mortality because of its antigenic variation. So far, very little is known about the antigenic patterns of the recent H3N2 
virus. By systematically mapping the antigenic relationships of H3N2 strains isolated since 2010, we discovered that two 
groups with obvious antigenic divergence, named SW13 (A/Switzerland/9715293/2013-like strains) and HK14 (A/Hong 
Kong/5738/2014-like strains), co-circulated during the 2014—2015 winter season. HK14 group co-circulated with SW13 in 
Europe and the United States during this season, while there were few strains of HK14 in mainland China, where SW13 has 
dominated since 2012. Furthermore, we found that substitutions near the receptor-binding site on hemagglutinin played an im- 
portant role in the antigenic variation of both the groups. These findings provide a comprehensive understanding of the recent 


antigenic evolution of H3N2 virus and will aid in the selection of vaccine strains. 
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The human influenza A (H3N2) virus has been one of the 
leading causes of seasonal epidemics since its emergence in 
1968 and poses a substantial threat to public health. H3N2 
infections appeared to become more severe during the 
2014-2015 flu season, which was widespread across the 
globe [1-4]. According to the Chinese Centre for Disease 
Control and Prevention (China CDC) [4], all of China expe- 
rienced H3N2 flu activity (Figure S1A). The percentage of 
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specimens positive for H3N2 was significantly higher during 
the 2014—2015 flu season than that during the previous four 
flu seasons (Figure S1B). 

As indicated in a report published by the US Centers for 
Disease Control and Prevention (CDC) [5], in 2014-2015 
the flu-associated hospitalization rate among people who 
were 65 years and older was the highest since the CDC be- 
gan surveillance in 2005. Another CDC report indicated that 
more than two thirds of the influenza A (H3N2) strains cir- 
culating in 2014-2015 were antigenically distinct from the 
H3N2 vaccine strain A/Texas/50/2012, which was used to 
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produce vaccines for the 2014-2015 season in the Northern 
Hemisphere [2]. Early estimates of this vaccine’s effective- 
ness showed the lowest efficacy of the last eight years (Fig- 
ure SIC). A deep understanding of the antigenic divergence 
of recent H3N2 strains is important for formulating proper 
surveillance and control strategies [6]. 

Recently, many efforts have been made to investigate the 
antigenic diversity, evolution, and geographical transmis- 
sion of the human influenza A (H3N2) virus [7—10]. It was 
reported that the antigenic evolution of the H3N2 virus was 
more punctuated than the genetic evolution [7]. A few ge- 
netic changes near the receptor-binding site (RBS) on he- 
magglutinin (HA) might significantly drive antigenic 
changes in the H3N2 virus [8]. Other researchers found that 
human influenza A (H3N2) strains originating in East and 
Southeast Asia usually seeded worldwide epidemics [9]. 
The most important factor driving the global dynamics 
of the human influenza A (H3N2) virus is air passenger 
flow [10]. Some very useful prediction models have been 
developed in recent years to accurately and rapidly catch the 
antigenic changes [11-14]. Developments include our own 
predict antigenic cluster (PREDAC) method [14], which has 
shown good performance in clustering the antigenically 
similar strains of the human influenza A (H3N2) virus. 

The PREDAC method provides a means to capture the 
antigenic diversity and explore the regional transmission of 
different groups of antigenically similar strains of the hu- 
man influenza A (H3N2) virus. By analyzing the antigenic 
relationships of all H3N2 strains isolated since 2010, we 
discovered that these strains could be divided into five anti- 
genic groups that showed the punctuated pattern of antigenic 
evolution. Unexpectedly, we found that two antigenic 
groups, SW13 and HK14, co-circulated during the 
2014-2015 flu season. 


1 Materials and methods 


1.1 Sequence data and phylogenetic analysis 


The HA1 protein sequences of 11,073 H3N2 viruses, sam- 
pled from all over the world between 2010 and 2014, were 
downloaded from the China CDC website and the Global 
Initiative on Sharing all Influenza Data (GISAID) platform. 
All HA1 sequences were aligned with Muscle v3.8.31 [21]. 
Re-assortment sequences from the laboratory and marked 
egg isolates were excluded. To reduce systematic biases, 
sequences with a gap content greater than 5% were re- 
moved. We excluded 25 sequences that were clearly recog- 
nized as outliers in the reconstructed strain trees. The phy- 
logenetic tree was constructed using PhyML [22] and dis- 
played in Dendroscope [23]. We collected structural data 
for H3N2 HA (A/Aichi/2/1968) [24] from the Protein Data 
Bank (PDB ID: 3HMG). 
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1.2 Mapping the antigenic groups of the human influ- 
enza A (H3N2) virus 


We used the computational method PREDAC [14] to model 
the antigenic groups of the H3N2 virus. The antigenic rela- 
tionship between each pair of strains was predicted using 
HA1 sequences, and pairs that were predicted to be anti- 
genically similar were grouped using the MCL(Markov 
Cluster) program [25]. The antigenic correlation network 
was visualized using Cytoscape [26], and the positions of 
some nodes were adjusted to distinguish the groups from 
each other. 


1.3 Molecular clock analysis 


The different phylogenetic relationships of the strains pre- 
sent in antigenic groups SW13 and HK14 were determined. 
We conducted a phylogenetic analysis for all strains within 
five years using the software BEAST v1.7.5 [18]. To reduce 
the number of very similar sequences in our initial data set, 
strains with the same sequence, country, and month of ob- 
servation were counted only once, and no more than five 
strains were selected from each continent in the same year; 
this reduced the number of strains to 145. Nucleotide se- 
quences encoding the HA1 subunit were downloaded from 
GISAID and China CDC. 

The Markov chain Monte Carlo analysis was performed 
in two independent runs. The codon position-specific HKY 
plus Gamma 14 nucleotide substitution model and the un- 
correlated lognormal molecular clock model were used in 
the analysis. The first 10% of trees were discarded as 
‘““burn-in’’, as recommended in the BEAST package manual. 
The tree with maximum log clade credibility was selected 
and visualized using FigTree v1.4.2 [27]. 


1.4 Percentages of strains in antigenic groups 


The percentage of strains in each group, on a monthly basis, 
was computed as the number of strains of one antigenic 
group divided by the total monthly number of strains. To 
investigate the antigenic patterns of different regions, the 
percentage was computed as the number of strains in one 
month divided by the number of strains of six months be- 
fore or after that month. 


1.5  Hemagglutination inhibition data and antigenic 
cartography 


HI data were collected from the European Centre for Dis- 
ease Prevention and Control [16] (see Table S1 for detailed 
information). We then used AntigenMap 3D [28] to create 
the antigenic map. 


1.6 Analysis of amino acid substitutions and N-glyco- 
sylation sites 


Cluster-difference amino acid substitutions were defined as 
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previously described [7], with a conservative cut-off of 90%. 
N-glycosylation sites were predicted using the NetNGlyc 
server (http://www.cbs.dtu.dk/services/NetNGlyc/). Bootstrap 
resampling (1,000 replicates) was performed. 


2 Results 


2.1 Antigenic evolution of the human influenza A 
(H3N2) virus since 2010 


The antigenic relationships of the human influenza A 


HK14 


Percentage of 
antigenic clusters 
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(H3N2) strains isolated over the last five years were ana- 
lyzed using the PREDAC method [14]. All of the seasonal 
human influenza A (H3N2) strains from 2010 to 2014 could 
be divided into five antigenic groups (Figure 1A): BRO7 
(A/Brisbane/10/2007-like strains), PEO9 (A/Perth/15/2009- 
like strains), VIC11 (A/Victoria/361/2011-like strains), 
SW13, and HK14. These groups were named for the vac- 
cine strain or representative strain they contained. 

Antigenic group BRO7 was replaced by PEO9, and PE09 
was replaced by VIC11. Two antigenic groups (SW13 and 
HK 14) co-circulated after the replacement of VIC11 in 


SW13 


3C.2 


HK14 | 3C.2a 


— VIC11 
— SW13 
— HK14 


Figure 1 Genetic and antigenic evolution of human influenza A (H3N2) virus during 2010-2014. A, Predicted antigenic correlation network and antigenic 
groups of human influenza A (H3N2) strains. The predicted antigenic groups are named according to the vaccine strain they contain. B, Phylogenetic tree of 
hemagglutinin subunit 1 (HA1) protein sequences of 2,221 strains. Subclades are colored and named as in (A). C, Antigenic cartography of isolates of anti- 
genic group HK 14 (belonging to clade 3C.2a) and antigenic group SW13 (belonging to clade 3C.3). Groups HK14 and SW13 are colored orange and red, 
respectively. D, Dynamic changes in the percentages of strains in antigenic groups on a yearly basis. The bottom bar represents the vaccine strains recom- 
mended by the WHO. The purple bar represents vaccine strain A/Brisbane/10/2007. The blue bar represents vaccine strain A/Perth/16/2009. The green bar 
represents vaccine strain A/Victoria/361/2011. Grey backgrounds indicate winter influenza seasons in the Northern Hemisphere. 
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March 2014 (Figure 1B and 1D). Based on a report from the 
World Health Organization (WHO) [15] and our phyloge- 
netic analysis, we concluded that both SW13 and HK14 
evolved from group VIC11 and fell into two distinct genetic 
clades, referred to as genetic clades 3C.2a and 3C.3 in the 
WHO report (Figure 1B). To confirm the antigenic distinc- 
tiveness of groups SW13 and VIC11, we collected hemag- 
glutination-inhibition (HI) test data from the surveillance 
reports of EU/EEA (the European Union and European 
Economic Area) countries (Table S1) and constructed a 
three-dimensional antigenic map (Figure 1C). Although the 
two groups were still not antigenically distinguishable from 
each other according to a WHO report in March 2015 [15], 
obvious antigenic divergence was observed between these 
two groups (Figure 1C) [16]. 

A significant difference in emergence and transmission 
was observed between antigenic groups SW13 and HK14. 
Group SW13 was predominant after October 2013, 
while group HK14 was noticed in 2014, according to the 
WHO [17]. Group HK14 first appeared in September 2013 
and co-circulated with group SW13 after April 2014, repre- 
senting over 20% of circulating strains. 

To infer the phylodynamic spread of the two co-circu- 
lating antigenic groups in time, we performed molecular 
clock analysis using the BEAST software package [18]. The 
inference was made for HA1 segments of strains within five 
years, as shown in Figure 2 (the detailed tree is shown in 
Figure S2). Groups SW13 and HK14 were distinct from each 
other and both derived from group VIC11. The date for the 
most recent common ancestor for group SW13 was August 
2012 (95% Bayesian credible interval May 2012—November 
2012), while that for group HK14 was September 2013 
(95% Bayesian credible interval May 2013—January 2014). 


2.2 Epidemic patterns of human influenza A (H3N2) 
virus among different regions 


For better understanding of the detailed epidemic patterns in 
different regions, we mapped the antigenic groups of human 
influenza A (H3N2) virus from 2010 to 2014 in three dif- 
ferent regions: mainland China, the United States, and Eu- 
rope (Figure 3). Overall, the epidemic patterns of H3N2 in 
the United States and Europe were more similar to each 
other than to the pattern in mainland China from 2010 to 
2014. The antigenic group PEO9 was replaced by group 
VIC11 in the 2011-2012 winter season in mainland China 
but was epidemic until the 2012-2013 winter season in the 
other two regions. The strains of antigenic group SW13 
were observed sporadically in the 2011-2012 winter season 
and increased quickly after November 2012. 

During the 2012-2013 winter season, the strains of anti- 
genic group SW13 were widespread in mainland China, 
Europe, and the United States, with a proportion almost 
equal to that of the predominant group, VIC11. Later, in the 
2013-2014 winter season, group SW13 was predominant in 
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Figure 2 Temporal reconstruction of the phylogenetic relationships of 
HA1 sequences. Antigenic groups are shown on the right. Internal nodes 
are reconstructed common ancestors; the 95% highest posterior density 
intervals on their dates are indicated by the orange bars. 


all three regions. The interesting thing is that another anti- 
genic group, HK14, also occurred during the 2013-2014 
winter season in the United States and Europe. It co-circu- 
lated with group SW13 during the 2014-2015 winter sea- 
son. However, in mainland China, SW14 was the dominant 
group during the 2013-2014 and 2014-2015 winter seasons 
and only a few strains of group HK14 were observed. 


2.3 Antigenic and genetic variation of the co-circulat- 
ing antigenic groups 


The observation of two unusual co-circulating antigenic 
groups drove us to further investigate the genetic and anti- 
genic variation between groups SW13 and HK14. We found 
that both HK14 and SW13 possessed several amino acid 
mutations compared with group VIC11 (Figure 4A). Most 
mutations were located in antigenic epitopes. The transition 
from VIC11 to SW13 was mainly associated with two sub- 
stitutions at positions 128 and 142. However, the transition 
from VIC11 to HK14 was mainly associated with three oth- 
er substitutions, at positions 144, 225, and 311. 

We then mapped the mutations onto the HA structure 
(Figure 4B and 4C). Four of the five substitutions (all ex- 
cept the substitution at position 128) were in or near the 
RBS, which is thought to play an important role in antigenic 
evolution [8]. Differences in glycosylation sites between the 
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Figure 3 Comparison of antigenic patterns of human influenza A (H3N2) strains in different regions. Dynamic changes in the percentages of strains in 
antigenic groups were recorded on a yearly basis in mainland China, the United States, and Europe. Different-colored curves represent different antigenic 
groups. Grey backgrounds indicate winter seasons. Vaccine strains recommended by the WHO for Northern and Southern Hemisphere influenza seasons are 


shown at the bottom. 


antigenic groups were also examined (Figure 4D). In com- 
parison with group VIC11, there were two changes (glyco- 
sylation sites 126 and 144) in group SW13 and two changes 
(glycosylation sites 144 and 158) in group HK14. These 
glycosylation sites are all located in antigenic epitopes A 
and B, which are thought to contribute significantly to anti- 
genic change [8]. 


3 Discussion 


Through systematic prediction and analysis of recent anti- 
genic relationships, we have developed a comprehensive 
picture of the antigenic evolution of human influenza A 
(H3N2) strains from the 2010-2011 to 2014-2015 winter 
seasons. Five antigenic groups were identified, and two of 
them, SW13 and HK14, co-circulated during the 2014-2015 
winter season. We also found that the epidemic patterns of 
antigenic groups SW13 and HK14 were quite different in 
different regions. During the 2014-2015 winter season, 
HK14 co-circulated with SW13 in the United States and 
Europe, while mainland China was still dominated by 
SW13 and only a few strains of HK14 were reported. 
Previous studies have indicated that mutations near the 


RBS or in the antigenic epitopes of HA are the major de- 
terminants of antigenic variation of influenza viruses [8,19]. 
Several mutations near the RBS or in the antigenic epitopes 
were found in antigenic groups SW13 and HK14. The 
T128A and R142G substitutions of group SW13 are located 
on epitopes B and A, respectively. Among the three muta- 
tions of group HK14, N225D is located in the RBS, N144S 
is quite near the RBS, and Q311H is located on epitope C. 
These mutations probably explain the antigenic variation of 
influenza A (H3N2) strains during the 2014-2015 winter 
season. 

Unlike with human influenza B virus [20], co-circulation 
of two antigenic groups rarely happens with human influ- 
enza A (H3N2) virus. However, during the 2014-2015 
winter season, co-circulation of antigenic groups SW13 and 
HK14 occurred in many regions. In the United States and 
Europe, SW13 and HK14 were present at similar scales 
during the most recent flu season. 

It is worth noting that the SW13 vaccine strain was not 
recommended until the 2015 Southern Hemisphere influen- 
za season, and the HK14 group was not noticed until re- 
cently. Thus, there would be a mismatch between the 
2014-2015 vaccines and the dominant human influenza A 
(H3N2) virus strains. Antigenic groups SW13 and HK14 
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Figure 4 Variation in amino acid substitutions and N-glycosylation sites between antigenic groups. A, Amino acid mutations in antigenic groups VIC11, 
SW13, and HK14. Five randomly selected strains from each group are shown. The top row shows the type of mutation site (pink sites, epitope A or B; blue 
sites, epitope C, D, or E; grey sites, other regions). The numbering of sites is shown at the bottom. B, Amino acid substitutions mapped on the HA1 protein 
structure. Differences between VIC11 and SW13 and between VIC11 are HK14 are colored red and orange, respectively. C, Substitutions near the recep- 
tor-binding site (blue) of HA. D, Predicted N-glycosylation sites in antigenic groups VIC11, SW13, and HK14. Error bars represent the standard deviation of 


bootstrap analysis with 1,000 replications. 


should both be considered during the selection of vaccine 
strains. Furthermore, knowledge of the antigenic patterns in 
different regions and the phylodynamic spread of each 
group, in time and space, will help improve vaccine strain 
selection. 
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